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Abstract 

In this contribution we present an intrinsic description of time-variant Port Hamil- 
tonian systems as they appear in modeling and control theory. This formulation 
is based on the splitting of the state bundle and the use of appropriate covariant 
derivatives, which guarantees that the structure of the equations is invariant with re- 
spect to time-variant coordinate transformations. In particular, we will interpret our 
covariant system representation in the context of control theoretic problems. Typi- 
cal examples are time-variant error systems related to trajectory tracking problems 
which allow for a Hamiltonian formulation. Furthermore we will analyze the concept 
of collocation and the balancing/interaction of power flows in an intrinsic fashion. 
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1 Introduction 



Hamiltonian systems are the object of analysis for a long period and they have 
been investigated from many different points of view and in many different sci- 
entific areas. In the last two decades, in mathematical physics especially field 
theoretic aspects of Hamiltonian systems without control input are of impor- 
tance, see [5,6,7]. In field theory the use of bundles to distinguish dependent 
and independent coordinates is commonly used and since time- variant lumped 
parameter systems can be seen as a special case of field theory with only one 
independent coordinate the use of bundles also applies to time-variant sys- 
tems where the fibration is accomplished with respect to the time-coordinate. 
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Beside field theory, also in classical mechanics, especially in the time-invariant 
setting, the geometric interpretation of the Hamiltonian picture is well estab- 
lished, see for example [1] for many details concerning this subject. 
From the control theoretic point of view, the class of Port Hamiltonian systems 
are a well-analyzed class, see for example [11,16] and references therein, where 
both the theoretical point of view and, of course, the physical applications play 
a prominent role. Roughly speaking, the main idea of many passivity based 
control approaches is to maintain the Hamiltonian structure of the system 
by feedback since this structure has some pleasing properties concerning the 
stability proof also in the nonlinear scenario. 

In the literature most of these approaches for the lumped parameter scenario 
concerning control theoretic aspects present system analysis, modeling and 
control for time- invariant systems, whereas the time-variant case is analyzed 
very rarely. We believe that the main difficulty in the time-variant scenario 
is the fact that the geometric picture of the equations changes considerably. 
In the time-invariant setting the role of the time is solely to be the curve pa- 
rameter, which is not true in the time-variant scenario. Contributions which 
treat the time-variant case especially with regard to control theoretic prob- 
lems are for example [3,4] where the authors consider, what they call gen- 
eralized Hamiltonian systems and canonical transformations which might be 
time-dependent. 

Two important applications where time-variant systems arise quite naturally 
based on a time-variant change of coordinates should be mentioned at this 
stage: Firstly, the introduction of displacement coordinates with respect to a 
system trajectory as it arises for instance when the analysis of the tracking 
error is the objective, see for example [4]. And secondly, in mechanics/robotics 
floating/ accelerated frames of reference are commonly used with respect to an 
inertial one. 

The main contribution of this paper are that: (i) an intrinsic definition of time- 
variant (Port) Hamiltonian systems is given based on a covariant derivative 
induced by a connection; (ii) this intrinsic description is analyzed in a dif- 
ferential geometric way; (iii) for the system class of time-variant (controlled) 
Hamiltonian mechanics, a covariant version of the power balance relation in- 
cluding collocation is developed; (iv) for the special case where (beside a possi- 
ble feed-forward) the connection can be expressed as an additive Hamiltonian 
the results of [3] are recovered. 

It is worth mentioning that in our opinion a time-variant (Port) Hamiltonian 
system has to be introduced using covariant derivatives, which differs signif- 
icantly from the definition in [3,4]. We identify 'covariant' with the fact that 
system properties do not depend on the chosen coordinate chart, i.e., we for- 
mulate systems in an intrinsic way. The key idea is the use of a connection 
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which induces a covariant derivative, see [5]. Partially, results in this paper 
have been presented preliminarily in [13,15]. 



2 The Time-Invariant Case 

This introductory section is a reminder of time-invariant Port Hamiltonian 
systems [8,11,16] including also the arising matching conditions when state 
transformations and affine input transformations are considered. It serves as a 
basis for the generalization to the time-variant case and is also used to intro- 
duce the differential geometric language which is then extensively exploited in 
the time- variant scenario. The notation is similar to the one in [5], where the 
interested reader can find much more details about this geometric machinery. 

To keep the formulas short and readable we will use tensor notation and es- 
pecially Einsteins convention on sums where we will not indicate the range of 
the used indices when they are clear from the context. We use the standard 
symbol (g) for the tensor product, d is the exterior derivative, J the natural 
contraction between tensor fields and o denotes the composition of maps. By 
d% are meant the partial derivatives with respect to coordinates with the in- 
dices Q. 

To study the time-invariant case of Port Hamiltonian systems in a geometric 
fashion we introduce the state manifold X equipped with coordinates (x a ), 
where a — 1, ... , dim(A') and we consider diffeomorphisms (in the sequel also 
called transition functions) of the type x = <p(x) where x denotes the states 
in the transformed coordinate system. Standard differential geometric con- 
structions, see [1,5,9,12] lead to the tangent bundle T(X) and the cotangent 
bundle T*(X), which possess the induced coordinates (x a ,x a ) and (x a ,x a ) 
with respect to the holonomic bases d a and dx a . Typical elements of T(X) 
(vector fields) and T*(X) (1-forms) read in local coordinates as w = x a (x)d a 
and oo = x a (x)dx a , respectively. To introduce in- and outputs we consider 
the vector bundle U — > X with the coordinates (x a ,u l ) for U and the base 
for the fibres where i = 1, . . . , dim(Wj-), where Ujr denotes the fibres of the 
input bundle (vector spaces) as well as the dual output vector bundle y — >■ X 
possessing the coordinates (x a ,yi) and the fibre base e\ Greek indices will 
correspond to the components of the coordinates of the state manifold and 
induced structures. Latin indices correspond to the components of the input 
and the output variables (fibres of the dual bundles U — > X and y — > X). Let 
us consider the maps J, R : T*(X) — > T(X) which are contravariant tensors 
that are given by the local coordinate expressions 

J = J aP d a ®dfi, R = R a/3 d a ® dp (1) 

with J a P,R a P E C°°{X) where J is skew-symmetric, i.e. J af} = -J^ a and R 
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is symmetric R a/3 = R^ a and positive-semidefinite. Furthermore we introduce 
the bundle map G : U — > 7~(X) which is a tensor that has the local coordinate 
expression G = G"e l ® d a with Gf G C°°(X). Having the maps J, R and G at 
our disposal a time-invariant Port Hamiltonian system (with dissipation), see 
[8,11,16] can be constructed as 

x = (J-R)\dH + G\u 
y = G*\dH 

where the function H G C°°{X) denotes the Hamiltonian and G* : T*(X) — > y 
the adjoint (dual) map of G. The local coordinate expression of (2) reads as 

x a = (j afi - R afS ) d H + Gfu* 

V J p ( 3 ) 

Vi = Gfd a H. 

We want to analyze structure preserving transformations for the system (2). 
To allow for affine input transformations we can replace the input bundle by an 
affine one Z — > X (with underlying vector bundle U — > X), for the geometric 
properties of affine bundles see for example [5] and references therein. The 
transition functions for the vector bundle and the affine bundle read as 

u = Mu, vJ = M^u i (4) 
u = Mu + g, u j = M{u i + g j (5) 

with M/,^ G C°°(X) where u denotes the transformed input coordinates 
and we restrict ourselves to regular transformations (i.e. M is invertible). 
The geometric representation of the system leads to the observation that the 
structure of (2) is preserved by a diffeomorphism of the type x = <p(x) together 
with (4). The case of an affine input bundle is more challenging since the 
preservation of the structure demands to solve a partial differential equation. 
See also [2] in this context, where the problem of general feedback equivalence 
of nonlinear systems to Port Hamiltonian systems is discussed and so called 
matching conditions appear. 

Lemma 1 Consider the system (2) together with the diffeomorphism x = 
(p(x) and (5). The structure of (2) is preserved if an only if we can find a 
solution H G C°°(X) of the partial differential equations 

( _ R^)d B H - (d a <p s G?Migi) = o. (6) 

Here J a/3 and R al3 are the components of the transformed tensors (1) with 
respect to x — (p(x). The inverse maps are denoted by x = (p(x) and M/Mj = 
5? where 5 is the Kronecker delta. 
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Remark 2 The partial differential equations (6) are written in the coordinates 
x but it is readily observed that it can be formulated in the original coordinates 
x, as well. 

The proof of this Lemma is a straightforward calculation in local coordinates. 
If in Lemma 1 a solution for H can be obtained, then the following Corollary 
is an immediate consequence. 

Corollary 3 Suppose (6) is met, then the system (3) in the new coordinates 
reads as 

-a = ( Jap _ R*P}dp(H _ ft) + Gfu 1 

m — Gfda(H — H) 
Wi th Gj = (d a ip^GfMfj o and H = H o (p. The output 

y s = Mj(y i -(d a H)G?). (7) 
and the Hamiltonian are transformed affine and in general u l yi ^ u?y-j is met. 

Example 4 Let us consider the shifting of a nonzero (but constant) equilib- 
rium point (x e , u e ) of the system (2). This leads to affine relations of the form 
x = x — x e and u = u — u e . For a different interpretation concerning Port 
Hamiltonian systems with nonzero equilibrium consider [10] whereas in [8] the 
construction of Lyapunov functions is treated in this case. 



3 The Time- Variant Case 

Time-variant systems arise when the time coordinate is treated as an addi- 
tional coordinate in contrast to the case where it is solely a curve parame- 
ter. This has the consequence for (Port) Hamiltonian systems that the maps 
J, R, G and the Hamiltonian H may depend on the time coordinate and that 
also time-dependent coordinate changes have to be taken into account. Con- 
cerning time-dependent changes of coordinates in mathematical physics, es- 
pecially in mechanics a covariant treatment of the equations requires to use 
frames of references formulated on bundles, see [5] . 

Remark 5 A special case of a time-variant coordinate transformation is the 
introduction of displacement coordinates with respect to a system trajectory of 
the (time-invariant) system (3). 

We will apply intrinsic concepts to model Hamiltonian systems. To motivate 
how these covariant concepts arise let us consider a time- variant system mod- 
eled on an extended state manifold given by the direct product B x X that 
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includes the time coordinate t° for B (the index will always correspond to 
the time coordinate in the sequel) as it is analyzed for instance in [3] and by 
a slight abuse of notation we have 

d Q x a {t) = J a P{t,x{t))d p H{t,x{t)) + G«{t,x(t))u\t) 
Vi (t) = G?(t,x(t))d a H(t,x(t)). 

The crucial point is now how to interpret the time derivative dox a (t) either 
using the structure of jet spaces (jet bundles) which are affine or using tangent 
structures. Let us consider a transformation of the form x a = Lp a [x^ , £°) with 
f° = t°. We have 

£ a = d <pH° + dptpW (8) 

and consequently the vector field do is mapped to the vector field do + doip a d a - 
This follows from i° = 1 together with (8) and the main reason for that 
fact is the observation that a trivial product bundle structure of the form 
B x X — > B is not preserved by time-variant transformations and therefore 
not the adequate choice for time-variant systems. To overcome this problem 
the system has to be formulated on a bundle £ — > B. The choice of the 
trivial one B x X — >■ B corresponds to a specific trivialization £ = B x X, i.e. 
to the choice of a certain reference frame. An intrinsic approach demands a 
formulation on the bundle £ — > B and the choice of a connection, which can 
be seen as a reference frame. 



3.1 Geometric Setting 

To obtain an appropriate geometric picture for the time-variant case we con- 
sider the bundle £ — )■ B, where we use coordinates (t°,x a ) for £ and obtain 
the following geometric structures. The tangent bundle T(£) with coordi- 
nates (t°,x a ,t°,x a ), the cotangent bundle T*(£) equipped with coordinates 
(t°, x a , to, x a ), as well as the vertical bundle V(£) with coordinates (t°, x a , x a ) 
and V*(£) with (t°, x a , x a ), where V(£) possess the induced bases d a and since 
V*(£) does not possess a canonical base without the choice of a connection it 
is denoted as dx a at this stage, see [5]. 

To be able to deal with time- derivatives of sections s : B — > £, i.e x = 
s(t) we introduce the first jet manifold J l [£) with the adapted coordinates 
(t°,x a ,Xg), see [5,12]. The coordinates Xq are often called derivative coor- 
dinates and the transition functions with respect to the bundle morphism 
x a = {p a (t, x) , £~° = t° (we do not consider time-reparametrization) read as 
x% = doip 51 + dpip^xo 3 . 

The key object for a covariant system representation will be a connection 
together with the associated covariant differential. The main philosophy be- 
hind a connection is the fact that T{£) possesses a canonical subbundle, 
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namely the vertical tangent bundle V(£) but there is no canonical horizon- 
tal complement such that T{£) = V{£) © %{£) holds, unless one specifies a 
horizont af^kubbundle and this is exactly what a connection does. By du- 
ality one derives analogously a decomposition of the cotangent bundle as 
T*(£) = V*(£) ®H*(£), where in this case %*(£) is canonically given and the 
connection defines V*(£). 



3.1.1 Connection 

Given a bundle £ — > B whose fibration induces the vertical tangent bundle 
V(£) it is the desire to obtain a splitting of the form T{£) = V(£) © %{£). 
This can be achieved by a connection V which in local coordinates can be 
represented as a tensor of the form 

r = dt° © {d + r«d a ) (9) 

with the connection coefficients Tq G C°°(£), see [5,12]. The transition func- 
tions for T when no time-reparametrization (i.e. i = t) takes place read as 

r* = do<p* + r%d a <p s (io) 

where Tq are the connection coefficients in the transformed coordinate system. 

This connection (9) can be used to define a covariant derivative and to split the 
tangent bundle T(£) — > £, i.e. T{£) = V{£) © %{£) as stated above, where 
%(£) denotes the horizontal subbundle. We have the following coordinate 
representations for the splitting for a typical element of T(£) and T*(£), 
respectively 

i°d + x a d a = i^w^ + (x a - i°T^ d a 
i dt° + x a dx a = (t + x a T^ dt° + x a tu% 
with the vector field Wq and the 1-form ujy 

w ^ = d + T«d a , u% = dx a - T^dt (12) 

such that wff and Uy qualify as bases for %{£) and V*(£ ), respectively. 

Remark 6 It is readily observed that a time-variant transformation x a = 
ip a (t, x) converts a trivial connection Tq — in one coordinate chart to a 
non-trivial one in the transformed coordinates Tq = dof a , see (10). 

1 It should be noted that the concept of vertical and horizontal parts of several 
objects will become important in the sequel and we will indicate this using V for 
vertical and T~L for horizontal. However, Ti should not be confused with the Hamil- 
tonian H. 
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3.1.2 Covariant Derivative 

Given the connection (9) a covariant differential relative to the connection 
T can be introduced as a map V r : J l {£) — > T*(B) <g> V(£) which reads in 
coordinates as 

V r = (a£ - T°) dt° <g> d a . (13) 

If the connection is trivial, i.e. Tq = then the covariant differential on the 
bundle £ — >■ B corresponds to the classical time derivative. Based on the 
covariant differential one can define the covariant derivative of a section s : 
B — > £ which follows as 

v r (s) = (d s a -r>«) dt° ® d a . 

Given the vector field d : B — > T(£>), the contraction 

d \v T (s) = (d s a -r«o S )d a 

is said to be the covariant derivative of s : B — > £ along <9 , see [5]. 
3.2 System Representation 

We define a time-variant (Port) Hamiltonian system on the bundle Z — > £ — > 
B together with a connection (9) such that the structure of the system is 
invariant with respect to bundle morphisms of £ — > B that explicitly depend 
on the time coordinate but we do not consider time-reparametrization (i.e., 
x = (p(t,x) and i = t). The input bundle U ->■ £, (t°,x a ,^) ->■ (t°,x a ) 
is a vector bundle, and Z — > £ is used in the affine case, analogous to the 
time-invariant case. 

Applying a time-variant transformation does not destroy the structure of a 
(Port) Hamiltonian system if it is defined in an intrinsic manner. But in con- 
trast to [3], where only the relative motion with respect to a frame should be 
expressed in a Port Hamiltonian framework, it is necessary in general to adopt 
the involved differential operators to preserve the structure of the system. This 
will be done by applying covariant derivatives. The choice of T corresponds 
to the selection of a frame of reference, i.e. starting in an inertial frame (with 
trivial connection) a time-variant transformation induces a non-trivial con- 
nection. An example how a non-trivial connection can arise by changing the 
frame of reference for a mechanical system is given in Section 4.2. 

Definition 7 Given a connection (9) then a Hamiltonian system on a bundle 
Z — )> £ — >■ B is given by 

c? J V r = (J - R)\ dH + G\u (14) 
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with J, R : V*(S) — > V(£) where J is skew symmetric and R is symmetric and 
positive semidefinite. Additionally we have the bundle map G : U — > V{£) and 
the Hamiltonian H G C°°(£). In local coordinates we obtain the expression 

x « - = - R aP ) dpH + Gy\ 

where 

v a HV = - R a ^ d p H + Gfu* 

is referred as the vertical part of the Hamiltonian vector field vuy 
with H, Gf, J a/3 , R aP E C°°(S). 

To show that in the time- variant scenario J and R are maps from the vertical 
cotangent bundle to the vertical tangent bundle J, R : V*(S) — > V(£) we have 
to observe that the total differential of the Hamiltonian dH = doHdt + 
dpHdx 13 can be decomposed according to (11) as 

dH = dpHufi + w£(H)dt° (16) 

with the 1-form and the vector field Wq from (12) from which we deduce 
the result, together with (13). 

Remark 8 Since the exterior derivative d can also be decomposed into a ver- 
tical and a horizontal one d = dy + d^, see [5,12] one can equivalently use 
d\>H instead of dH in the relation (14)- 

Let us consider the total time change do(H) of the Hamiltonian H along 
solutions. For a time-invariant system we obtain a relation of the form 

do(H) = v H (H) = -(d a H)R a ?(dpH) + u'yi 

where vh{H) denotes the Lie-derivative of the Hamiltonian H along the Hamil- 
tonian vector field vh = — i? a/3 ) dpH + G"u\ see (3) and the total time 
derivative (time-invariant case) reads as do = x a d a . 

The intrinsic version in the time- variant case is given in the following corollary. 

Corollary 9 The total time change of the Hamiltonian leads to a decomposi- 
tion of the form 

d (H)=w H (H) + v H , v (H) (17) 

or in local coordinates 

do(H) + T«d a {H)-{d a H)R a P{d p H) + d a (H)G?u l 

S v * v ' 

w^(H) v H ,v(H) 

2 It should be noted that = dx a — T^dt serves as an adapted basis for V*(£) 
induced by the connection T. This justifies J, R : V*{£) — > V{£) as we have claimed. 



(15) 
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where the vector field w H is induced by the connection T, see (12), the vector 
field vh,v from (15) and do = do + x^d a . It is worth noting that we have a de- 
composition of the derivative operator into a horizontal vSq and a vertical vh,v 
component where the horizontal part degenerates to d only if the connection 
is trivial, i.e. Tq = 0. 

This follows by a direct calculation based on the intrinsic system representa- 
tion. The discussion of collocation will be performed using a special class of 
systems in the next section. 

3.3 A special class of Port Hamiltonian Systems 

Mechanical systems in a Port Hamiltonian representation are distinguished 
since the bundle £ — > B has an even richer geometric structure. The main 
difference is the separation of the x a coordinates in positions q a and momenta 
q a . To obtain the correct geometric picture we introduce the bundle Q — > B 
with coordinates (t°,q a ) for Q. This bundle structure induces again some 
tangent structures where we make use of the dual vertical cotangent bundle 
V*(Q) — > Q with coordinates (t°, q a , q a ). Furthermore we will utilize T(V*(<2)) 
which possesses the adapted bases (do, d a , d a ) with d = and to be conform 
with most of the literature we make the identification q a = p a . 

3.3.1 The Composite Bundle Structure 

From the bundles introduced so far V*(<2) — > Q and Q — > B we can construct 
the composite bundle structure V*(<2) — > <2 — B. 

The additional fibration V*(Q) ->■ B, (t°,q a ,p a ) ->• (t°) plays the role of the 
state bundle £ — > B. We choose a connection 7 that corresponds to the selec- 
tion of a frame of reference that splits T(<2) 



Based on the connection 7 the so-called covertical connection, see [5] 



can be constructed that splits the tangent bundle T(V*(<2)) with respect to 
the fibration V*(<2) — > B and this is the connection that corresponds to T, see 
(9) that splits T{£). 

In classical mechanics it is well known that the Hamiltonian vector field can be 
defined using Symplectic and Poisson structures, see [1,5]. These concepts can 
also be generalized to the time-variant scenario. We base the construction of 



7 = dt° ® (do + 7o a £U , To e C°°(Q). 



(18) 




(19) 
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the Hamiltonian vector field using the Hamilton form uh which follows from 
the canonical Liouville 1-form (by a pull-back) 

u H =p a dq a -(H + Pal «)dt° 

with H e C°°(V*(Q)). 

From the relations vh\ duin = and vh\ dt° = 1 the autonomous Hamiltonian 
vector field v H : V*(Q) ->■ T(V*(Q) follows with # 7 = H + as 

Uh = d + d Q # 7 d a - <9 Q # 7 d Q . (20) 

Remark 10 If the connection 7 is trivial, i.e. 7o = holds then the Hamil- 
tonian vector field reads as vh = d + d a Hd a — d a Hd a since then — H 
which is the standard result in mechanics, see [1]. 



3.3.2 Power Balance Equation for Controlled Mechanical Systems 

The connection (19) enables us to split the Hamiltonian vector field (20) into 
a vertical and horizontal part, respectively, see also [13]. It follows that (anal- 
ogously to (11)) we have a decomposition of the form 

v HV = d a Hd a -d a Hd a 

(21) 

vh,h = do + To d « - {d P lQ)vpd p 

according to the Hamiltonian connection (19). Let us inspect the relation (17) 
in the case of a mechanical system where it is obvious that now the vector 
field Wq from (12) is replaced by vh,h due to our richer bundle structure. 
We consider the autonomous case first and since no dissipation is present we 
obtain 



v H (H)=v H , H (H) (22) 

where vh(H) denotes the Lie-derivative of the Hamiltonian H with respect 
to the vector field vh- To include control inputs we consider an extended 
Hamiltonian corresponding to a controlled Hamiltonian, see [9], of the form 

H = H - H c y (23) 

with H , H C)P e C°°(V*(Q)) and the input functions up e C°°(B). 

Based on (23) and the relation (22) we can state the following Theorem con- 
cerning a covariant formulation of the power flows as well as the concept of 
collocation. 
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Theorem 11 The change of the free Hamiltonian H along solutions of the 
Hamiltonian system decomposes as 

v H (H ) = v HtH (H ) + v H:V (H CtP )u p 

where we used the controlled Hamiltonian H = H — H c p u p and the decompo- 
sition (21). 

Obviously the choice of the output y p = vh,v{H C}P ) allows a physical inter- 
pretation of the power flows of the system, since vh,h(H ) corresponds to the 
power caused by the free Hamiltonian H and the product y p u p describes the 
power flow into the system caused by the input. Theorem 11 is an intrinsic 
version of the balancing/interaction of power flows and it reduces to the well 
known formula when the connection is trivial, i.e. 7g = 0. It is worth mention- 
ing that the splitting of the Hamiltonian field as in (21) is essential to obtain 
a coordinate free representation. 



4 Applications 

This section is devoted to a short discussion of two possible applications of 
the presented theory. Firstly, we focus on the description of the error system 
arising when the stabilization of the trajectory tracking error is the objective 
in a Port Hamiltonian framework and discuss the role of the connection in 
this context. As a second application we present the equations of motion of 
a mass particle observed from a rotating frame of reference, in a time-variant 
Hamiltonian formulation since this example demonstrates how to calculate 
the non-trivial connections coefficients and additionally one can show that 
this connection is used to formulate conservations laws and/or power balance 
relations. 

4-1 Trajectory tracking 

This section deals with time-variant Port Hamiltonian systems as they arise 
quite naturally when a feed-forward based approach is applied to a time- 
independent Port Hamiltonian system and the control objective is the stabi- 
lization of the error system using techniques of passivity. It is evident that the 
error system is constructed using a time- variant transformation which enables 
us to discuss the developed machinery on this concrete example. It should 
be pointed out that compared to the approach in [4] we use a different defi- 
nition of a time-variant Port Hamiltonian system (see Definition 7) which is 
coordinate independent and more general due to the connection term. If the 
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connection is interpreted as an additive term in a modified Hamiltonian then 
our results coincide with the results of [4]. 

Let us consider a Port Hamiltonian system of the form 

The control system (24) is modeled on a bundle S — >■ B where the coordinates 
are obviously adapted to the connection, such that the connection coefficients 
then read as Tq = 0. If a desired trajectory Cd(t) and a corresponding input 
T]d(t), which produces this trajectory, is given 

d c a d = ((J a/3 - R aP ) d p H + Gf Vd ) o Cd (25) 
then a transformation with respect to a reference trajectory can be stated as 



x a = ^(t , a J J ) = ^-cg) (26) 
u^MUu'-vi), M}eC°°{X). (27) 

Remark 12 The relations (26) and (27) are of course only a trivial choice 
leading to an error system. But for our purposes, the geometric interpretation 
of time-variant Port Hamiltonian systems, the exact structure of (26) and (27) 
is of minor importance, since the intrinsic system formulation is independent 
of the special choice of the transition functions. For a more general discussion 
of error systems in the Port Hamiltonian context, see [4]. 

It is readily observed that the connection coefficients read as 

f * = dotf = -{docX. (28) 

with ip a from (26). 

The combination of the results of Lemma 1 together with the fact that in the 
time- variant case a connection appears that leads to an additional affine term, 
we obtain the following Corollary. 

Corollary 13 Given Q(t) and rj d (t) fulfilling (25) then the transformations 
(26) and (27) applied to the system (24) lead to a representation as 

-a = gap _ R«P)dp(H + H)+ Gfu\ (29) 

if and only if the partial differential equations 

(J^-R^)dpH= f^ + gG^ 

allow a solution for H. 
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From Corollary 13 we deduce that if the connection * and the feed-forward part 
** are expressed as an additive Hamiltonian, one obtains a Port-Hamiltonian 
representation which is beneficial when the control objective is to stabilize 
error systems that arise typically when trajectory tracking is the demand 



4-2 Rotating Frame of Reference 



To show that a time-variant transformation preserves the Hamiltonian struc- 
ture we will consider the equations of motion of a mass particle using two 
different frames, i.e. an inertial one and a rotating one with respect to the 
inertial frame. Let us consider an inertial system with Euclidean coordinates 
(g Q ,t°) together with the bundle structure Q — > B. The canonical equations 
of motion for a mass particle with mass m G K + read as 



d s a = d a H , d (p a o s ) = -d a H (30) 
with s : B — > Q. For this example the Hamiltonian is given as 

where 6 denotes the Kronecker delta. 

A rotating coordinate chart with respect to the inertial one can be constructed, 
using (f = Rl(t°)qP = <p s (q p ,lP) with Rp a/3 R^ = 5^. The non-trivial con- 
nection in the floating reference system can be computed as in (10) and reads 
as 7? = d (R%)R§q? = Qfq?. 

The equations of motion in the rotating coordinate system follow as, see (20) 
or [14] for a detailed exposition concerning covariant derivatives 

/0 _ (31) 

d (Pa s) +p / 5<9a7o = -d a H . 

Comparing (30) with (31) it is evident that they both are Hamiltonian repre- 
sentations of the same physical problem, but in contrast to the inertial frame 
where the left hand side consists of partial time derivatives only in the rotating 
frame a differential operator induced by the connection (19) has to be applied. 
In the concrete example where the connection reads as 7o = ^q' 5 together 
with s : B — ?■ Q the equations follow as 



d s s -^=-p^ (32) 
3)(p a os)+j^ = 0, (33) 
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when no other forces are applied. 

Remark 14 To apply Theorem 11 it should be noted that vg n has to be 
constructed with the connection 7q = . Then the Hamiltonian H is a 
conserved quantity as the conservation law vg n (H) = (see Theorem 11 with 
Hq = H ) is met. If we apply forces to control the system then the full relation 
of Theorem 11 has to be applied including the additional expression vg v , see 
(21). 

Combining the equations (32) and (33) we obtain 

d (d ^ - np?)5^ + (d tp - ny)5^ni = o 

which can be rewritten as 

doos* - d npP - 2ftjflbs* + = 

where we have exploited the skew symmetry of the so-called angular velocity 
tensor Q. This is a classical result in mechanics, i.e. this is an expression for 
the acceleration including the Coriolis and the Centrifugal acceleration. The 
relations (32) and (33) are the covariant version in a Hamiltonian point of 
view. 

Summarizing we can state, that time-variant transformations preserve the 
Hamiltonian structure, if a (Port) Hamiltonian system is introduced covari- 
antly as in Definition 7. The system properties can be expressed in a time- 
variant frame which may lead to non trivial connections, which then requires 
the use of covariant derivatives. Also the conservation laws and the power 
balance laws have to be formulated in an intrinsic fashion, see Theorem 11. 

Finally it should be stressed again, that the interpretation of the connection 
in the subsections (4.1) and (4.2) is completely different. In subsection (4.1) 
the connection is absorbed in a modified Hamiltonian, whereas in subsection 
(4.2) the connection is explicitly part of the covariant differential. 



5 Conclusion 

In this paper we investigated the geometry of time-variant (Port) Hamilto- 
nian systems and used an intrinsic description based on connections and co- 
variant derivatives. We were interested in two concrete applications where 
time-variant systems arise quite naturally, namely error systems with regard 
to trajectory tracking in control theory as well as time- variant mechanics. Con- 
cerning control theory we described that the formulation of the error system 
in a Hamiltonian fashion requires to absorb the nontrivial connection and the 
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feed-forward part in a modified Hamiltonian. Here the key problem was the 
definition of a time-variant Hamiltonian system itself, i.e. if the tensors are 
time-dependent and/or the coordinate chart is moving. The interpretation of 
these two scenarios is significant for a correct understanding of time-variant 
systems. Furthermore time-variant mechanics has been analyzed in a covari- 
ant way, where we showed how collocation and the balancing/interaction of 
power flows can be formulated in an intrinsic way. 
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